Dipolar confinement-induced resonances of ultracold gases in waveguides 
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Exploring the impact of anisotropic dipolar interactions for ultracold collisions in quasi-one- 
dimensional traps we derive and analyze dipolar confinement-induced resonances (DCIRs) attributed 
to different angular momentum states. Our analysis is based on an extended A-matrix approach for 
quasi-one dimensional geometries. An analytically derived resonance condition allows us to predict 
the positions of DCIRs. The results obtained within the A'-matrix theory are in excellent agreement 
with corresponding exact numerical calculations. 
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In low-dimensional geometries due to tightly confin- 
ing traps, ultracold atomic scattering undergoes crucial 
modifications yielding the effect of confinement-induced 
resonances (CIRs) [3, A variety of confinement- 

induced processes have been demonstrated recently 0- 
e.g. fermionic p-wave CIRs @, dual-CIR @, res- 



or excited channels and 
Hoj . The deepened un- 



onant molecule formation 7[ 
higher partial wave CIRs 9, 
derstanding of CIR physics has lead to major achieve- 
ments in the experimental manipulation of in- 

ter acting gaseous atomic matter. Together with the 
extensive study of free-space dipolar collisions flil - 1231 ]. 
confinement-induced resonant scattering introduces an 
intriguing perspective in the control of dipolar many- 
body phases [24[. Indeed, reduced dimensionality has 
lead to the control of the internal and external degrees 
of freedom of molecules allowin g th e realization of dense 
ultracold polar molecule gases [25[. In view of the sub- 
stantial theoretical effort made on confined dipolar scat- 
tering 2^- 31 1, the need for a rigorous understanding of 
the role of anisotropic forces in CIRs becomes evident. 

In this letter, we analytically derive resonance con- 
dition for s-wave dipolar CIRs (DCIRs), with an ex- 
plicit dependence on the dipole-dipole interaction (DDI) 
strength. This is done within an extended AT-matrix 
formalism for harmonic quasi-one-dimensional (Q1D) ge- 
ometries [l(| , which incorporates the DDI and takes into 
account contributions from higher angular momentum 
states. These £-wave states are firstly coupled due to 
the anisotropic nature of the DDI and secondly by the 
harmonic confinement, which leads to a rich resonance 
structure of the DCIRs. Interestingly, this interplay be- 
tween the confinement and the DDI leads to an intricate 
dependence of the s-wave DCIRs positions on the dipo- 
lar interaction strength. This is in contrast to the usual 
case of s-wave CIRs [1| , where their positions is given by 
a s /a± = 0.68 with a s being the s-wave scattering length 
and <xj_ the length of the transversal confinement. The 



derived resonance condition thus reveals in detail the im- 
pact of the DDI anisotropy on the CIR effect and provides 
the necessary tool for the experimental control of dipolar 
collisions in Q1D. 

Setting the scene, we consider a system of two bosonic, 
nonreactive polar molecules which collide in a Q1D 
waveguide. They are treated as perfect dipoles fully 
polarized by an external electric field along the waveg- 
uide axis z. The transversal confinement is induced 
by a two-dimensional (2D) harmonic potential yielding 
a separation of the center of mass and relative degrees 
of freedom. The physics of the collisional processes is 
then completely described by the relative Hamiltonian 
H = — |^V 2 + ^w^p 2 +V r i n t(r) expressed in cylindrical co- 
ordinates r = (p, <j), z), where fi denotes the reduced mass 
and wj_ is the confinement frequency. The interaction po- 
tential is 14 n t(r) = V sh (r) + - 3(z ■ r) 2 ], where the 
short-range term V s h is modeled by a Lennard- Jones (LJ) 
12-6 potential, 14h(r) = ttt ~ 7^; an d the second term 
describes the DDI where d is the induced dipole moment. 
The justification of the LJ-potential for the short-range 
behavior of the two nonreactive molecules lies in the sep- 
aration of the energy scales for elastic ultracold colli- 
sions and inelastic chemical processes [3l|. The ranges 
of the LJ and DDI potentials are Z v dw = (2uC6/ft 2 ) 3 
and Id = iid 2 /H 2 , respectively. As in Refs.[B|,[I(|, we con- 
sider the confining oscillator length a± = y/h/fiu!± to be 
the larger length scale, i.e. a± ^> l v dw. Id- 

This condition separates the configuration space into 
three domains with respect to the relative distance r: 

(I) Id, Ivdw < T <C a±: At small separation distances 
the interactions Mnt(r) dominate, so that the dipoles ef- 
fectively experience a free-space collision of total energy 
E = h 2 k 2 /2p, with the according symmetry imposed. 
In this region, the corresponding wave function can thus 
be efficiently expanded in the £-wave angular momentum 
eigenstates. Note that, due to the azimuthial rotational 
symmetry of H, the quantum number m is conserved, 
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and will here be fixed to m = 0. On the contrary, the an- 
gular momentum £ is not conserved due to the anisotropy 
of the DDI, which couples states with A£ = £ — £' = 2. 
The scattering information of Vint is then imprinted in a 
free-space A-matrix of tridiagonal and symmetric form 
in ^-representation, with entries for even £ due to bosonic 
symmetry; considering up to <?-wave contributions, it 
reads 





(K ss 


K sd 





A 3D = | 


K ds 


K dd 


Kdg 




V o 


K gd 


K gg 



Since the waveguide confinement is not experienced in 
this region, the entries of A 3D are a measure of the dis- 
tortion of the wave function only by the presence of Vi n t . 

(II) a± <C r — > oo: At large separations the cylindrical 
confinement prevails, such that the wave function decom- 
poses into the radial eigenmodcs n of the transversal 2D 
harmonic oscillator. These are regarded as asymptotic 
channels of the scattering process, which is described by 
the corresponding A- matrix in one dimension, A 1D . 

(III) l d , l v dw <C r <C a±: At intermediate separa- 
tion distances there is, in general, an admixture of £- 
wave states from Vint and n-mode states from the con- 
finement. However, in this particular domain, both the 
interaction and confining potentials essentially vanish, so 
that the corresponding scattering solutions can be effi- 
ciently matched. Therefore, a non-orthogonal local frame 
transformation U_ exists [32j, which enables the projec- 
tion of the wave function of domain (I) onto that of do- 
main (II). In particular, U_ transforms the corresponding 
A-matriccs into one another: K = U 1 K OJJ U with 
its elements being (K_ 1D ) nn > = J2w U^KwUtn' . Note 
that both A 3£> and K_ 1D depend on the total energy E. 

In the following, we consider low-energy scattering in 
the asymptotic transversal ground state n = 0, follow- 
ing the criterion qoa± <C 1 for the relative longitudinal 
wave vector go which is considered to be small and fi- 
nite. The total collision energy, in domain II, is E = 
H 2 k 2 /2fi = hbj± + h 2 q 2 /2/j,, so that only the first channel 
(n = 0) is energetically open (o), while higher channels 
(n > 0) remain closed (c), since E < haj±(2n + 1). The 
asymptotic components of the wave function in the c- 
channels, however, contains exponential divergences, ren- 
dering the scattering unphysical. This behavior is reme- 
died in the framework of multichannel quantum defect 
theory (MQDT) by imposing the physically acceptable 
boundary conditions in c-channels [33[ , yielding a physi- 
cal A"-matrix given by 

A oo = Koo + % Koc (I - iKcc ) Kco , (2) 

where I is the identity matrix. The roots of det(X — 
iK_l®) = provide the bound states in the c-channels 



which energetically lie in the continuum of the o- 
channels. 

Using the Dyson- like form of K_ OQ from Ref. we 
obtain 



— 00 " q a^{l + a.-ny K ' 

where the scalar products are short-hand for summed 
coefficients of the a,£ and 7L 1 given explicitly in the 
supplementary material (SM). £ and contain traces 
J2 n U^gUi'n over the c-channels n, yielding combinations 
of the Hurwitz C-function for each (£, ^')-pair. ol con- 
sists of single-term combinations of a U ' = a u > /a± , where 
a ee > = —K u > jk are generalized, energy dependent scat- 
tering lengths for dipolar collisions in free-space [20| ■ Due 
to Eq. (fT]) we remark that a contains up to g- wave a~w- 
terms. Eq. © therefore encapsulates directly the impact 
of the anisotropy of the DDI on the confined scattering: 
A 00 is determined by the dipolc-induced a u < -terms con- 
tained in a, which arc simultaneously weighted by the 
coupling due to the confinement via £ and TZ. 

The resonant behavior arises in the form of poles of 
K_ oa , given by the roots of the equation 

1 + 0.-71 = 0. (4) 

These coincide with the zeros of det(Z — iK l J?) = in 
Eq. ©, which demonstrates that the origin of the res- 
onance structure is a Fano-Fcshbach mechanism. Note 
that the resonance condition, Eq. ([?]), can be met in mul- 
tiple ways by allowing either one of the aw -terms to be 
dominant. 

It is thus evident that due to the presence of the DDI 
within the waveguide, different £-wave states from the 
domain (I) contribute to this resonance mechanism. The 
£-wave labeling is still used in the sense that a partic- 
ular partial wave dominates over the others, although 
they are coupled together due to DDI, which leads to 
broad s-wave and narrow (£ > 0)-wavc DCIRs with corre- 
sponding positions in the parameter space, e.g. l d , l v dw, 
a±, determined by Eq. Note that this difference of 
the widths arises from the free-space dipolar scattering, 
where within the adiabatic approximation of the two- 
body dynamics [22[ the s-wave adiabatic channel pos- 
sesses an effective — 1/r 4 potential tail, while each £ > 
state exhibits repulsive barriers leading to increasingly 
narrower resonances with increasing €-wave character. 

We will now focus on the s-wave DCIRs, which 
are of immediate experimental relevance. Thus, we 
isolate the free-space a ss dipolar scattering length 
on the left hand side of Eq. ([!]), and we obtain 
the corresponding resonance condition a ss (ka± , d) = 
J"({a«'(/ca_L, d)}, {TZi(ka±)}), where 
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T = - 



1 + addj^ + Qgg^-12) + agg(R-3 + a s dR-ii) + flsd(7?4 + ad g TZi 3 ) + adgCR-5 + adgT^-i) + a^ d (TZ 6 + aggTZu) 

Hi + d dd Kw + a gg Tl$ + a dg TZ 9 + (a dg ~ d dd a gg )Hu 



(5) 



with the 7?-i (i = 1 . . . 14) being the components of 7Z. 
(see SM). Originating from the poles of K_ oa in Eq. ([3]), 
the quantity J- now embodies the interplay between 
the DDI anisotropy and the confinement. Note that in 
general a ss (ka±,d) and J~({dw(ka±, d)}, {lZi(ka±)}) de- 
pend differently on the dipole moment d, thusly their 
equality provides a transcendental equation being ful- 
filled for particular values of d in the parameter space. 

The s-wave DCIRs appear in the vicinity of the free 
space resonances which arc s-wave dominated, i.e. near 
the broad divergences of the a ss dipolar scattering length, 
where all the higher aw scattering lengths are non 



resonant. Additionally, as was shown in Ref.[2CJ, [27 1 
within the Born approximation (BA) the non-resonant 
aw, except £ = i 1 = terms are proportional to the 
dipolar length l d . Consequently, this universal, i.e. V s h- 
independent, behavior of the higher aw terms leads to 
the following simplification of Eq. ([5]) : 



i + vik + mil + mi 3 a 

CO + dJd + a4\ 



(6) 



where l d = l d /a± and the r]j and Oj consist of combi- 
nations of TtiS (see SM). In the considered low-energy 
limit qoa± <C 1, they acquire the values r\\ ps 1.844, 
rj 2 w -1.119, ?73 ~ 0.013, a « -1.46, crj « 2.008 
and <72 ~ 0.046. For l d = 0, the resonance condition, 
a ss = Tba, reduces to a s = — l/oo = 0.68, as expected 
for the s-wave CIR [l[. 

The analytical results from the BA, Eq. are com- 
pared in Fig[T]to those of Eq. {[5J> , in which the aw are cal- 
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FIG. 1. (color online) T versus l d = l d /a±. The black line 
refers to Tba (Eq. (|6j)), whereas J> (red line) and T s (blue 
line) refer to Eq. (J5]) where the terms a U i are numerically 
calculated for a LJ potential that possesses for d = a weakly 
or a strongly s-wave bound state, respectively. 
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FIG. 2. (color online) The transmission T: analytical re- 
sults (solid line) and numerical calculations for (a) a 3 3> Z v dw 
(— O — ), (b) a s <C /vdw ( — □—). (c) and (d) show the corre- 
sponding quantity \det(X — iK^ c )|, for the parameter values 
of (a) and (b), respectively. 



culated numerically from the dipolar free-space problem 
34|. We see that, whereas the BA breaks down close to 

free-space resonances (aw — > oo, £, £' > 0) as expected 
35|, J-ba is in good agreement with the non-resonant 

parts of T\ these are also the regimes of interest here, 

since the s-wave DCIRs occur away from aw, ^ £' > 

free-space resonances. 



To study the universal aspects of J 7 , we consider for 
d = the two limiting cases of a weakly and a strongly s- 
wave bound state in the LJ potential, yielding a s 3> l v <rw 
and a s -C ^ v dw, respectively. We observe in Fig. Q] that 
T s (strongly bound-blue line) and T w (weakly bound- 
red line) practically coincide everywhere apart from the 
positions of the narrow resonant features. The non- 
universality of these resonant features mainly arises from 
the coupling term a s d, which is strongly affected by 
the non- universal and strongly resonant behavior of a ss . 
However, the distance between the corresponding diver- 
gences in F s and .F w decreases as T d increases. This is 
because, as the DDI becomes stronger, it dominates the 
short-range LJ interaction and eventually shields it com- 
pletely, thereby restoring the universal behavior of a ss . 

We now investigate the resonant structure of the trans- 
mission coefficient T for confined dipolar scattering, 
which we derive analytically in terms of the physical K- 
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FIG. 3. (color online) (a) a ss depicted for a s S> l v dw (— O — ) 
and a s <C lvdw (—□ — ), whereas the black line shows the 
quantity Tba, (b) the analytically calculated transmission T 
for a s 2> ^vdw (red line) and a s <C lvdw (blue line), (c) and 
(d) represent magnifications of the gray shaded areas of (a) 
and (b), respectively. 



matrix (see SM) as 
T 



my 



(7) 



where the numerically calculated a. from the unconfincd 
problem is used as an input in Eq. (|3J), as before. The 
transmission is shown in Fig(5]as a function of Id, again 
for the cases of (a) a weakly and (b) a strongly bound 
s-wave state in the LJ potential at d = 0, each featuring 
sequences of s- and higher partial waves (£ > 0) DCIRs. 
As mentioned previously, the latter are much narrower 
due to the presence of the repulsive barrier. These results 
are compared to exact numerical calculations of T (red 
circles and blue squares in Figj2] (a) and (b) , respectively) 
based on a scheme presented in Ref. [36[. An excellent 
agreement is observed. Moreover, FigJ5] (c) and (d) show 
the corresponding expression \det(I — iK^)\, which is 
observed to tend to zero exactly at the positions of T « 0. 
This indeed illustrates that the £-w&ve DCIRs fulfill a 
Fano-Feshbach scenario. 

FigJ3] (a) shows a graphical solution of the resonance 
condition for the positions of the s-wave DCIRs: the 
dipolar scattering within the waveguide becomes reso- 
nant for the values of Id where a ss = J- ~ J-ba > 0. We 
focus on the broad divergences of a ss related to s-wave 
dominated dipolar free space resonances, whose width in- 
creases with Id, since the DDI becomes more attractive. 
The intersections of Tba and a ss in Fig(3] (a) are seen to 
occur exactly at the zeros T(ld) = of the transmission in 
Figj3](b), as clearly demonstrated for the magnified grey 
shaded regions in FigJ3] (c) and (d). This striking coin- 
cidence demonstrates the high accuracy of the derived 



resonance condition for s-wave DCIRs, for both limiting 
cases a s <C l v <iw and a s > l v dw- 

This accurate prediction of the s-wave DCIRs is re- 
tained over the whole range of the DDI strengths, from 
the non-universal to the universal regime (Id > 0.11) of 
the s-wave dominated dipolar free space resonances. In 
the universal regime the short-range physics represented 
by V s h indeed becomes irrelevant, in the sense that the 
positions of the s-wave DCIRs (intersections of a ss and 
•7"ba(^)) for different strengths of V s h to a good approx- 
imation coincide in the (a ss , y-planc. 

FigJ3] (a) comprises the most intriguing difference of 
DCIRs compared to well-established CIRs with isotropic 
interactions: In contrast to the case of CIRs, the posi- 
tions of s-wave DCIRs, measured in values of a ss , in- 
creases for increasing interaction strength of the DDI 
potential, or, equivalently, for an increasing number of 
s-wave dominated dipolar free space resonances which 
have become bound within V5 n t(r). Therefore, in the 
presence of DDI the confinement-induced shift of free- 
space resonances is no longer constant, but increases for 
successive resonances passing the open channel threshold. 
This behavior arises from the anisotropic nature of the 
DDI, which mixes higher partial €-waves more strongly 
for larger Id, and is enhanced via their recoupling by the 
confinement, yielding substantial contributions to the po- 
sitions of the s-wave DCIRs. 

In conclusion, we have extended the if-matrix formal- 
ism to treat dipolar collisions and include higher angular 
momentum partial waves in the presence of a harmonic 
Q1D confinement leading to the prediction of £-wave 
DCIRs. In particular, we analyzed in detail the case of s- 
wave DCIRs and derived analytically the corresponding 
resonance condition in the form a ss — J-ba, containing 
explicitly the dependence on the DDI strength. Apart 
from providing an essential ingredient for the resonant 
control of nonreactive polar molecule gases in Q1D, this 
result sheds light on the physics underlying the DCIR 
effect: We demonstrate how the DDI couples the £-wave 
states of the short-range interaction potential, which are 
in turn mixed by the confinement. As a result of this 
interplay between DDI and confinement, the positions 
of the s-wave DCIRs, measured in values of a ss , were 
shown to shift for increasing strength of the interaction 
potential, in contrast to the case of mere s-wave CIRs. 
The theoretical advance presented here, combined with 
MQDT theory [37l | . may pave the way for new insights 
on reactive polar molecule collisions in quasi-2D confining 
traps and add to the understanding of the corresponding 
experimental realizations pjj . Furthermore, we remark 
that the s-wave dominated dipolar free space resonances 
can be induced either by strong dc-electric or laser fields 

H 
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This together with the tunability of the con- 
finement frequency provides us with excellent tools for 
probing s-wave DCIRs. 
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In the following we present the explicit form of the physical K 1 -matrix and the cor- 
responding scalar products, as well as the relations of the ^-constants which refer to the 
resonance condition and finally, we discuss the derivation of the transmission coefficient T. 



I. THE EXPLICIT FORM OF THE PHYSICAL K 1D -MATKIX FOR QUASI-1D 
GEOMETRIES 



Let us provide some insights into the derivation of the physical K} D '-matrix, namely , 
for dipole-dipole collisions. In general we have 



£Z = £2 + &z - izzr'ig, (i) 

where we assumed that only one energetically open (o) channel, namely n = and n ^ 
energetically closed (o) channels. Here K l ® = J2ee' ^>o' U refers to the open-open channel 
transitions and the matrix element K l ® ,u is defined as {K l ® ,u ) 00 = (U T ) eK u >Ui' with 
being the corresponding matrix element of K 3D . = Eee'lCc'"' (iCo = T,wlCo' U ') 
refers to open-closed (closed-open) transitions and each Klf ,ee (K x ® ,u ) is a single row- 
column-) matrix, where its entries are defined as (K_l® ,u )on — (U T )oeKu>Ue> n ( {K^' u ) n0 = 
(U T ) n eK U rU£'o ). K l ® = J2ep — cc' U re f ers to closed-closed transitions, and the matrix 
elements of K^' U ' are given by {K^ ,U ') n n' = (U T ) n iK ie iUi> n >where each Kj£' u ' is a rank 
one matrix. 

The matrix (I — iK 1 ®) = X — i J2w ^IF'" * s inverted analytically where for each K x ® ,u 
rank one matrix we use recursively the inversion identity for rank one matrices of Ref. 
Thus takes a Dyson-like form by rewriting Eq. © in terms of K^' u ' , K^ M \ K^' u ' 
and the analytically inverted X — i Y^et' — cc* ' U ^ OT eac h ££'— pair. Thereafter, performing 
analytically all the corresponding matrix multiplications Eq. ([T]) then reads 

-id = a • g 

-°° q a x (l + cx-ny {) 
where the scalar product cz • 1Z has the following form: 

a • 1Z = a ss TZi + a dd 1Z 2 + a gg TZ 3 + d sd H± + a dg K 5 + a 2 sd 1Z e + a 2 dg R 7 + a ss a gg TZs 
+ (i ss d dg Tlg + a ss d dd Rio + d gg d sd Rn + addfl g gTli2 + (i sd d dg Tli3 

+ (ttsstfig + tigged ~ ttss&dd(lgg)'TtlAi (3) 
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with a u > = a u > / a_i_ with a u > = —K u > / k being generalized, energy dependent scattering 
lengths for dipolar collisions in free-space [2j encapsulating the scattering information of the 
short-range and the dipole-dipole interactions. The coefficients TZi, i = 1, . . . , 14, are the 
components of the vector 1Z and are defined by the following: 
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where 6j with i = 1, . . . ,6 are constants given in the following: 
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here £(•, •) denotes the Hurwitz zeta function and e = 1/2 — (kaj_/2) 2 . 
Similarly, the scalar product ck ■ £ results in the following expression: 
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+ (a ss a^_ + (Xgg(X s( i OjssO'ddOjgg)^, 



14, 



(6) 



with £j, i = 1, . . . , 14, being the components of the vector £ and are defined as follows: 
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with the constants b{, i — 1, . . . , 6, given in Eq. (j5]). 

Note that in the low-energy regime, qoa± <C 1, the dimensionless total energy becomes 



II. THE 77- AND cj-CONSTANTS OF THE RESONANCE CONDITION 

The explicit form of the 77- and a-constants expressed in terms of the components of 1Z 
take on the following appearance: 
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III. THE EFFECTIVE ID HAMILTONIAN AND THE TRANSMISSION COEFFI- 
CIENT T 



The physical K_ 1D matrix and the corresponding scattering wave function in the asymp- 
totic regime, |r| — > oo, are related as follows: 



^ s (z,p) = [cosfobW) -lL D sin(g k|)]$o,o(p) (9) 



where \I/ phys is the physical wave function |3| and $0,0 (p) is the eigenfunction of the two- 
dimensional harmonic oscillator for n = m = 0, i.e. the lowest transversal channel. Note 
that \|> phys consists only of the real part of the total wave function in the asymptotic regime. 
Thus in order to relate the transmission coefficient T to the matrix we map the relative 
Hamiltonian H of the two dipoles in the waveguide onto an effective one-dimensional (ID) 
Hamiltonian where they interact via a ID delta function, V(z) = gioS(z). This is permitted 
since the relative Hamiltonian H asymptotically is fully separable, and all the relevant 
scattering information is imprinted in the ^-component of the wave function. Consequently, 
by imposing the boundary condition of the 5 function we determine the prefactor such 
that the real part of the wave function of the effective ID Hamiltonian is the same as the z- 
dependent part of \[/ phys . Thus the effective coupling constant g±o becomes gm = —^-^-K cp 
and the corresponding transmission coefficient T of the effective ID Hamiltonian reads 

T= 1 ~T^ = ^T7^- ( 10 ) 
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